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Abstract 

We study the effect of parameters uncertainties on a stochastic diffusion model, in 
particular the impact on the pricing of contingent claims, thanks to Dirichlet Forms 
methods. We apply recent techniques, developed by Bouleau, to hedging procedures 
in order to compute the sensitivities of SDE trajectories with respect to parameter 
perturbations. We show that this model can reproduce a bid-ask spread. We also prove 
that, if the stochastic differential equation admits a closed form representation, also 
the sensitivities have closed form representations. 

We exhibit the case of log-normal diffusion and we show that this framework foresees 
a smiled implied volatility surface coherent with historical data. 

Keywords: Uncertainty, Stochastic Differential Equations, Dynamic Hedging, Er- 
ror Theory using Dirichlet Forms, Bias, Bid- Ask Spread. 
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1 Introduction 

The purpose of this article is to study the sensitivity of the flow of a stochastic differential 
equation with respect to its parameters and, in particular, the adjustments on the computa- 
tion of conditional expectation. The main application is to evaluate the correction of option 
pricing due to the uncertainty on parameters driving the SDE that describe the prices. 

The usual framework of asset pricing is due to Black and Scholes [6], it provides condi- 
tions on the market which guarantee the existence of a hedging strategy to cover a contingent 
claim and, therefore, the existence and the uniqueness of its price. However, the hypothe- 
ses of classical asset pricing have clearly some deficiencies and drawback: they assume in 
particular a complete knowledge of the market and, in particular, of the parameters driving 
the SDE that describes the assets prices. 

The risk assessment is often carry out in terms of sensitivity with respect to a determin- 
istic variation of model parameters. These sensitivities are usually called greeks. However, 
this approach encounters mathematical difficulties, in dealing with infinite dimensions, and 
practical ones. A large literature exists in the subject to cover the risk of uncertainty. The 
most important example is the uncertain volatility model, see Avelaneda et al. [2], Lyons 
|26] and Barles [5J, which takes into account the difficulties to calibrate the volatility in 
the Black Schoels model. Their strategy is to apply the stochastic control to super-hedge 
a contingent claim under model uncertainty. We also recall the recent paper of Denis and 
Martini |14| where a more general framework is introduced to study super-replication and 
model uncertainty, their methodology make use of the Choquet capacities. 
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The principal drawback of the super-replication techniques is that the super-hedging 
cost is too high and the corresponding super-hedging strategy is too conservative. As an 
example, Kramkov |25| shows that the super-hedging cost corresponds to the value of the 
option under the least favorable martingale measure. Likewise, Bellamy and Jeanblanc [5] 
prove that in a jumps-diffusion model the price range for a call option corresponds to the 
interval given by the no-arbitrage conditions, e.g. the super-hedging strategy of a call option 
is to buy one unit of the underlying. It seems therefore that super-hedging is not an effective 
methodology, since the option buyer will prefer to buy the underlying rather than pay the 
same price to have an option. The option writer has to take some risk in order to propose 
a competitive price. In other words, the option writer has to take the risk to lose money, 
with a small probability that this event occurs, in order to have the opportunity to make 
money with the contract, i.e. selling the option and hedging it. 

Many attempts to attack the problem are appeared in literature, Avelaneda et al. [2] 
bound the volatility between two levels. Other ways are the pricing via utility maximization, 
see for instance El Karoui and Rouge |17| . or via local risk minimization, see for instance 
Follmer and Schweizer |18j . 

An alternative way, based on potential theory (for reference see Albeverio [1], Bouleau 
and Hirsch [7] and Fukushima et al. |20j), has been suggested by Bouleau [8], this method 
is an extension of Malliavin calculus. We assume that the uncertainties are supposed to be 
infinitesimal random variables in accord with statistical estimation. These random variables 
are generally poorly known, so we assume that only the two first order are known, i.e. the 
biases and the variances/covariances, since the accuracy is usually computed using statistical 
measures, like Fisher information. If the uncertainty on parameter is small, we may neglect 
orders higher than the second, i.e. we take into account only the bias and the variance of 
the parameters. This cut is justified by the evidence that the parameters estimation is often 
hard, owing to the non-linearity, the complexity of the equations and to the bad quality or 
the shortage of the data. In error theory using Dirichlet forms, we associate the variance to 
the "carre du champ" operator T and the shift to the generator of the semigroup A. The 
two operators are closed in a suitable space and have a closed chain rule, see Bouleau [9]. 

The main result of this paper gives closed forms for the corrections on the price of 
contingent claims due to the uncertainties on the parameters. These are corrections for 
an investor who searches to compensate the systematic bias and that accepts a residual 
risk on the parameters sensitivity. We define this risk taking in a statistical way. A direct 
consequence of this result is the separation of the buy and sell prices due to the asymmetry 
between the buyer and the seller with respect to the risk on parameters. An outstanding 
result is that a systematic bias exists even if the parameters are unbiased, this bias is due to 
the non-linearity of the payoff. Our analysis covers a large class of stochastic process with 
continuous path. In particular, equity models with local or stochastic volatility. 

We exhibit, at the end, an example in the case of Black Scholes model with uncertainty 
on volatility parameter. We show that that our analysis foresees a smiled implied volatility 
and We give the explicit formula for vanilla prices and the bid-ask spread. 

The paper is organized as follows: In section 2, we present a survey of error theory using 
Dirichlet forms technique. Section 3 is devoted to the study of the impact of uncertainty on 
a diffusion model. We analyze the profit and loss process and we compute its law depending 
both on the underlying diffusion and the parameter uncertainty. We also introduce a pricing 
principle to over-hedge the contingent claim and we exhibit the bid and ask prices. In section 
4, we give an example with log-normal diffusion without drift. We exhibit the bid and the 
ask prices and we prove, under some hypotheses, that the implied volatility has a smiled 
behavior. 
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2 Mathematical Framework for Uncertainty 



We begin by giving a general introduction to the study of the sensitivity with respect to a 
stochastic perturbation and a formal definition of the framework that we will use in accord 
with the error theory using Dirichlet forms in accord with Bouleau [8|. In this survey, we 
follow [9]. 

We consider a function F{U\, U2, ■■■), depending on parameters (Ui, U2, ••■) that we 
suppose afflicted with uncertainties. We assume that the function F is regular enough and 
we search to evaluate the impact on F of the uncertainties on U. The first study of this 
problem is due to Gauss that prove the following expansion for the variance of F if the 
uncertainties are small compared with the parameters values: 



However, this relation is proved only if the number of the parameters is fixed and if F has 
an explicit formula, i.e. if it does not be defined via a limit. The main aim of the error 
theory using Dirichlet forms is to go beyond these limits in order to evaluate the effect of a 
stochastic fluctuation on more complex objects like stochastic integrals. For that, we study 
a parameter x with a small uncertainty y/eY, on which we compute a non-linear function /. 
Our parameter x is replaced by the random variable X = x + \J~eY with variance e Var[Y]. 
We now apply the Taylor expansion to function / and we find 

Bias[/(X)] = E[f(X) - f(x)} = f{x) Bias[X] + ^f'(x) Var[X] + o(e) 

Var[/(X)] = E [{/(X) - f(x)} 2 ] = [f'(x)] 2 Var[X] + o(e). 

If we suppose e really small we can cut the high terms on e and find two closed chain rules 
for the bias and the variance. These are known in literature since the variance follows 
the same rule of a carre du champ operator of a probability space equipped with a local 
Dirichlet form, while the bias verifies the rule of the generator of the semigroup associated 
to the Dirichlet form, see for instance Bouleau and Hirsch [7j. The main advantage of this 
comparison is that the carre du champ and the generator of semigroup are closed operator 
with respect to the graph norm, see for instance Fukushima et al. |20j. Therefore, they are 
good operators to study objects defined by limits, like stochastic integrals. 

The axiomatization of this idea is introduced by Bouleau [8] as follows: He defines an 
error structure as a probability space equipped with a local Dirichlet form owning a carre 
du champ. 

Definition 2.1 (Error structure) 

An error structure is a term ^f2, F, P, B, T^j , where 



dF dF 



F, Pj is a probability space; 
1 is a dense sub-vector space of I? ( Q, F, P); 



r is a positive symmetric bilinear application from Dxl into L 1 \ Q, F, PJ satisfying 

the functional calculus of class C 1 Pi Lip, i.e. if F and G are of class C 1 and Lips chitzian, 
u and v £ B 7 we have F(u) and G(v) £ B and 

(1) T [F(u), G(v)} = F'(u) G'{v) T[u, v] P a.s.; 

The bilinear form £[u, v] = [T[u, v]] is closed; 
We generally write T[u] for T[u, u]. 
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With this definition £ is a Dirichlet form and V is the associated carre du champ oper- 
ator. The Hille Yosida theorem, see for instance Albeverio |T] and Fukushima et al. |20j . 
guarantees that it exists a semigroup, and, then, a generator A coherent with the Dirichlet 
form £. This generator (A, T>A) is a self-adjoint operator that satisfies, for F G C 2 , u G T>A 
and Y[u] G L 2 (P): 

(2) A [F{u)] = F'(u) A[u] + -F"(u) T[u\ P a.s.; 

moreover, it is a closed operator with respect to the graph norm. We underline two impor- 
tant results related to this theory. First of all the concept of error structure is deeply related 
to statistics. As a matter of fact, the uncertainties on parameters come from a statistical 
estimation and Bouleau and Chorro |10| have proved a connection between error theory 
using Dirichlet forms and Fischer information theory. Furthermore, the error structures 
have nice properties, in particular is possible to prove that the product of two or countably 
many error structures is an error structure, see Bouleau 

The main drawback of the carre du champ operator is its bi-linearity that makes com- 
putations awkward to perform. An easy way to overcome it is to introduce a new operator, 
the sharp. We recall the definition of sharp operator associated with T, see Bouleau and 
Hirsch [7], section II. 6. for the proof. 

Proposition 2.1 (Sharp operator) 

Let (^l, F, P, ID, be an error structure and ^0, F, a copy of the probability space 

(^l, F, P^ . We assume that the space D is separable, then there exists a sharp operator ( 
with these three properties: 



Mu L 2 (l 



• G D, T[u] = I (u*) 

. V« G D" and F G C 1 n Lip, (F(u\, ... , u n ))* = £™ =1 (§§- 



,# 



o uj u t 

The sharp operator is an useful tool to compute T because the sharp is linear whereas 
the carre du champ is bilinear. In analogy with the classical approach of error theory 
we associate the carre du champ operator T to the normalized variance of the error, the 
sharp operator becomes a linear version of the standard deviation of the error. Similarly, 
the generator describes the error biases after normalization, for more details we refers to 
Bouleau [8] chapters III and V and [9]. 

In order to apply the result of error theory to our problem, we assume that we can define 
an error structure for each parameter of our model, we denote them aj, and that we suppose 
the following assumption held for each parameter. 

Assumption H-l (Error theory framework) 

1. ai G T>A, T[ai] and A[ai] are known; 

2. the functions x i— y T[ai](x) and x i— > A[a>i](x) belong to L 1 n L 2 are continuous at 
x = ai and not vanishing}]; 

3. the error structure admits a sharp operator denoted ( )#. 



1 We use an abuse of language, since when we write T[a,i] and .4[a;], we mean that ai is a random variable. 
Afterwards, we compute F[a;] and A[a,i] on the point x = ai, where ai denotes the estimated value of the 
parameter. 
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We conclude with a direct corollary of this theory. The error theory using Dirichlet forms 
restricts its analysis to the study of two first order errors propagation, i.e. the bias and the 
variance. This fact is justified by the knowledge on the parameters uncertainties, generally 
given by the Fischer information matrix, that is often bad. The study of the high orders is a 
very hard problem for both mathematical and practical reasons. From mathematical point 
of view, we have to study the chain rules of the high orders, like skewness and kurtosis, and to 
prove that the related operators are closed in a suitable space. However, the crucial problem 
remains to have fine estimations for the high orders of the uncertainties, this obstacle cannot 
be overcome easily. Therefore, we decide to restrict our study to the two first order. As a 
consequence, we implicitly work with gaussian random variables. Finally we can state: 

Remark 2.1 (impact of uncertainty) 

The impact of uncertainty on the parameter transforms a constant into a gaussian dis- 
tribution of the form 



(3) F(X) ~ F(x) + e A[F{X)](X = x) + y/l y/F[F(X)](X = x) G 

where G is a standard Gaussian variable and e is a vanishing parameter. 

In a more conservative way, i.e. if we suppose that the gaussian approximation is too 
rough, we can apply the Chebyshev's inequality that state in our case: 

Proposition 2.2 (Chebyshev's inequality) 

The random variable F{X) verifies the following inequality for all k > 1 and supposing 
that e is a vanishing parameter. 



(4) P F(X) - F(x) - A[F(X)](X = x) > k y/eF[F(X)](X 



< 



1 



1 + k 2 



These results explain the role of the generator and the carre du champ operator, see 
Bouleau for a more general analysis. The theoretical image is perturbed due to the 
uncertainty on the parameter, this effect is small, however it produces not only a noise but 
also it alters the mean. 



3 Diffusion model under perturbation 

We start with the classical Black Scholes model, afterward denoted BS, see Black and 
Scholes [B]. Let (0, J 7 , P) be the historical probability space and Bt the associated Brownian 
motion, the dynamic of the risky asset under historical probability P is given by the following 
diffusion in accord with the model of Black and Scholes: 



(5) dS t = Stfidt + S t cr dB t 

In this framework, the price of a European vanilla option is well known. This model presents 
many advantages, in particular the pricing depends only on volatility and we find closed 
forms for premium and greeks of vanilla options. Unluckily, the BS model cannot reproduce 
the market price of call options for all strikes at the same volatility, this effect is called 
smile. To take into account this phenomenon, we analyze two main extensions, the local 
and the stochastic volatility models, afterward denoted LV and SV, see for instance Dupire 
|16] . Hull and White |23| and Heston [22] , In this two classes of models, the parameter ctq 
is replaced by a function a that depends on the time t, on the underlying St and, in the 
case of SV models, on a random source. The stochastic differential equation verified by the 
price of the underlying is 

(6) dS t = S t fidt + S t a(t, S t , u>) dB t 
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This class of model is large enough to include a lot of stochastic model currently used in 
finance, like SABR and fast mean reverting SV, see Hagan et al )21| and Fouque et al |19] , 
However, all diffusion with jumps model are excluded from our analysis, see for instance 
Cont and Tankov |13] , 

Henceforth, we denote (f2, J 7 , IP) the probability space where the Brownian motion is de- 
fined, Tt the standard filtration generated by the Brownian motion Bt and E the expectation 
under the probability P. 

Our goal is to analyze the sensitivities of these models with respect to the fluctuations 
on their parameters. As a matter of fact, all models depend on some parameters, generally 
few, that depend on the underlying and enable the calibration of the model. When we 
sort out a given value for a parameter, thanks to a calibration methodology, it is known 
with an uncertainty. This uncertainty can be estimated using statistical methods, like 
Fischer information, or computing the sensitivity of the model calibration. We underline 
that each method, both statistical and calibrative, yields parameters with uncertainties 
and these uncertainties have a random nature. Our goal is to analyze their impact on the 
management and the hedging of a contingent claim on the side of the seller that tries to 
escape any risk. 

We propose to consider a perturbation of the model given by SDE ([6]) by means of an 
error structure on the volatility function a(t, St, oj). 
We make the following hypotheses: 

Assumption H-2 (Asset evolution and uncertainty impact) 

1. the real market follows the SDE |6|) with fixed but unknown parameters, i.e. the un- 
derlying follows the SDE and this diffusion does not suffer the uncertainty on the 
parameters. The market is viable and complete, i.e. they are enough traded assets to 
guarantee that any contingent claims admits an hedging portfolio; 

2. the option seller knows that the underlying follows the SDE |6j) but does not know the 
function a, i.e. the values of the parameters; 

3. the option seller has to estimate the parameters of his model, so its volatility con- 
tains intrinsic inaccuracies, we model this ambiguity by means of an error structure; 
nonetheless we assume that the stock price St is not erroneous, but the price and the 
greeks of the option yes. The option seller evaluates the impact of this uncertainty 
on his profit and loss process, afterward denoted P&L, and tries to compensate it 
modifying his prices. 

Assumption H-3 (Uncertainty on volatility) 

1. the uncertainty on the volatility function a(t, St, w) is described by an error structure 
(f2, J-, P, D, r) associated to the function a; 

2. the uncertainty on a is independent of the remain of the model, i.e. of the probability 
space (Q, J 7 , P) where the Brownian motion Bt is defined; 

3. the error structure and the function a verify the assumption \H- 11 

4- the function a (t, x, oj) is square integrable and admits a series expansion, i.e. a(t, x, lo) = 
Si a i^i(i) x i w )> where the functions (pi belong to C > on (t, x) and are square in- 
tegrable on the third variable. We assume that the uncertainty is carried by the co- 
efficients ai, i.e. we set the aj to be random and the error structure is given by the 
product of the error structures on each a%. We also suppose held the property J^.2. in 
Bouleau f8$ chapter V page 84, thanks to that only a finite number of coefficients is 
non zero. 



6 



Before continuing, we discuss briefly the two previous assumptions. We have divided 
them depending on the nature of the assumption, more financial for the first one and 
mathematical for the second. Assumption IH-31 sets the mathematical frame used in this 
paper, points 2 and 4 are assumed to make easier the proof of results and to simplify some 
computation, but they can be replaced with more general hypothesis. 

The main financial hypotheses are resumed in assumption IH-21 The first point states 
that the underlying follows an exact SDE without perturbation. The uncertainty is added 
by the option seller when he searches to know the parameters of the SDE. This assumption 
is easy to understand in financial frame and avoids many problems on the mathematical 
framework. In particular when the volatility is uncertain the set of probabilities describing 
the whole class of possible probabilistic views is not dominated. In this case, the classical 
approach used in mathematical finance cannot be followed, the main examples about this 
subject are Denis and Martini |14] and Denis and Kervarec |15j . 

We also add the following hypothesis on the class of contingent claim analyzed in this 
work. 

Assumption H-4 (Contingent Claims) 

Let $ be the payoff of the contingent claim. We assume that <3? depends only on St and 
belongs to C 2 . We also suppose that the two first derivatives are bounded. 

It is plain that the previous assumption is very restrictive, since call and put does not 
verifies for instance. However, it is easy to define a series of payoff $j that dominates the 
payoff of a call (or a put) and converges uniformly and to take the limit using the hypothesis 
of no free lunch with vanishing risk. 

3.1 Management of profit and loss process 

We consider that the option seller uses his data and proprietary view to set the model 
parameters using an optimization procedure, see for instance Dupire |16j and Cont Tankov 
|13] , Given these parameters, the risk neutral- measure Q exists and is unique, this proba- 
bility enable him to define the fair price of the option and the hedging strategy under the 
hypothesis that we knowledge the real values of the paprameters. 

We study the profit and loss process associated to this hedging portfolio. We neglect, in 
this first work, that the option holder can sell the option, so we assume that he holds the 
contingent claim until the maturity. We also assume that all prices are denominated using 
the risk-free asset as numeraire. The profit and loss process at the maturity of the option 
seller is given by the price of the option plus the hedging strategy minus the final payoff 
that the seller has to pay to the holder, so we have 

(7) PkL{T) = F(q, x, 0) + [ T A( ? , St, t)dS t - $(5 T ), 

J o 

where F and A denote respectively the price of the option and its first derivative with 
respect to the underlying. The key remark about the P&L equation (J7|) is that the price 
and the delta depend on the volatility process, estimated by the trader, <j that is not the 
real market volatility a, i.e. the parameters used by the option seller to compute the price 
and the hedging position are their estimated parameters that are different from the real 
parameters in diffusion ([6]). It is plain that a good calibration produces a set of parameters 
very close to the real values but the uncertainty remains, as an example the variance of a 
unbiased statistical estimator is at least high as the Cramer Rao bound. Moreover, thanks to 
a result of Bouleau and Chorro |10| . V is equal the inverse of the Fisher information matrix, 
so the use of error theory using Dirichlet forms is efficient since it exploits all information 
on data given the fact that the Cramer Rao bound is verified. 
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We concentrate our attention on the law of the P&L at the maturity. In absence of 
uncertainty the random value P&iL(T) is equal to zero almost surely, i.e. the option can be 
exactly hedged. However, the exact hedging strategy cannot be performed since the option 
seller does not know the parameters value. Therefore, P&zL(T) is a random variable and 
we have the following remark. 

Remark 3.1 (Random sources) 

The value of profit and loss process at the maturity is a random variable depending on 
two random sources: 

• First of all, the stochastic "real" model J 7 , P) since the trader cannot use the 
correct hedging portfolio. 

• Second, the space (fl, J 7 , ¥), i.e. the stochastic process that depends on a random 
component independent to the brownian motion Bt- 

Remark 3.2 (Role of historical probability) 

The profit and loss process must be studied on historical probability P. As a matter of 
fact, the risk neutral probability Q can be used if and only if the market is complete, i.e. if 
all contingent claims are attainable. In our case, the option seller does not know the real 
diffusion of the underlying, so the law of the P&zL is not degenerate. The main impact is 
that the drift of the SDE f6jj plays a crucial role and that the second term in P&iL process 
is not a martingale. This fact complicates the computations on this paper. However, the 
important role of the drift fj, in asset pricing when the market is incomplete is emphasized 
in literature, let us mention Avelaneda et al. [2], Karatzas et al. 1241 and Lyons [26]. 

The two previous remarks show that the computation of the law of the P&lL process is 
an hard problem. Moreover, an essential requisite lacks to define the law of the PSzL since 
the PSzL depends on ? and the law of this process is poorly known since it depends on a 
calibration methodology. We assume, implicitly on assumption IH-31 that the option seller 
can estimate the mean value of his parameters and their variance, i.e. the two first orders 
of the law of <j. The knowledge of the high orders is a statistical problem often too hard to 
solve. 

Given the knowledge of the two first orders of the law of we can estimate the two first 
orders of the law of P&iL at best. This remark justifies our recourse to the error theory 
using Dirichlet forms, see Bouleau [8]. More precisely, we have remarked that the profit 
and loss process depends on two random sources, one is the stochastic process ? depends 
on some parameters that are random variables on the space (O, J-, P), the second is the 
Brownian motion Bt that live on the probability space (Q, J 7 , P). We can define accurately 
the law of PSzL with respect to the random source on the probability space P but we can 
just estimate the two first orders of the dependency with respect to P. 

In order to analyze the law of P&L process, it is sufficient to study the P-expectation 
on a class of regular test functions h(P&zL) and the error on them using error theory. In 
practice we will compute the bias and the variance of E[/i(P&L)], where E denotes the 
expectation under the probability P. To simplify our results, we suppose moreover the 
following assumption held. 

Assumption H-5 (Expansion Approach) 

We assume that the volatility estimated by the option seller is unbiased, i.e. we compute 
the value of our expectation for <j = a. 

The previous assumption involves in particular that the P&L is worth zero P-almost 
surely, so this hypothesis enable us to simplify our results in two ways. First of all, we don't 
need to distinguish the two volatilities a and we are interested only on the stochastic 
correction, in this sense we have an expansion approach. Secondly under P the random 
variable P&iL(T) is a constant. Then we can state the following theorem. 
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Theorem 3.1 (Law of the profit and loss process) 

Under assumptions \H-^ \H-3l H~4 and \H-5i we have the following bias and variance: 



(8) 
(9) 

where 



A\E[h(PkL)]] 

r [e [h (p&l)}} 



h'(0)A 1 (a) + -h"(0) A 2 (a) 



\ti 



9(a), 



dF 1 



d 2 F 



da(<pi 



;(q,x, 0) 2 (O, x) T[ai](a,i 



A 2 (a) 



dA 

a 2 A 



.A[oi](oi) 



;(s,S t , t) $(t, S t , L))S t 



dt r[oj](oi) 



T 



E 



OA 



\da{<t>i 



■(<?, 5s, s) ^(s, S s , w) ) .S 2 a 2 



ds 



*(<7) 



+ E 



\dcr((f>: 

dF 

da(4>i] 



0) 0,(0, cc) + /U 



(ct,x,0) <fo(0, x) + /i 



<9A 



■(<?, 5 S , s)<pi{s, S s , u)S s ds 



E 



dA 



(cr, S s ,s) S s (j>i(s, S s , lo) 



ds\ T[ai](a,i 



where g^.) denotes the Gateaux derivative with respect to a variation of the volatility along 
the component (pi . Moreover, we have the following truncated expansion 



(10) E [h(P&L)\ ~ h(0) + e ti(0) Ai(cr) + e ^"(0) A 2 (a) + [/i'(0)] 2 *(^) G, 

where G is a standard gaussian random variable. Otherwise, we have the following Cheby- 
shev's inequality 
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E [h(P&L)] - /i(0) - e ^ [E [h(PkL)]] > k y/e T [E [h (P&L)]] 



< 



1 



1 + yfc 2 ' 



Proof: We start with the study of the variance. Thanks to assumption IH-ll included 
in assumption IH-31 the operator T admits a sharp operator, in particular in each sub-error 
structure associated to each parameter dj it exists a sharp operator denoted ()#. Moreover, 
the error structure, defined in assumption IH-3| is the product of the sub-error structures on 
each parameter, so the different sub-error structures are independent among themselves. 

We now study the sharp of E \h (P&iL)]. Thanks to the linearity of the sharp operator 
and the expectation and the smoothness of the test function h we have 



(E[h(PkL)]) # = E 



(12) 



E 



h'{P&zL) (P&L)* 
h! {P&L) |f( ? , x, 0) + J' A( ? , S t , t)dS t - §{S T ) 



#■ 



Using the linearity of the sharp on braces, we remark that (<&(5r))# = since the payoff is 
independent on the volatility estimated by the option seller. The same argument and the 
closedness of the sharp enable us to write 



£ A( ? , S t , t) dS t ^j = £ (A( ? , S t , t)f dS t 
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where we have used the first point of assumption IH-21 Now, we recall that the price and the 
hedging position on the underlying are smooth functions with respect to the volatility since 
the SDE ([6]) has a gaussian kernel and thanks to assumption IH-41 Therefore, we can take 
the Gateaux-derivatives of F and A with respect to the volatility change in the direction of 
the functions fa, denoted and we have 



(13) 
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dA 
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■(?, S t , t) fa{t, S t , u)dS t 



where we have make use of the expansion of the volatility given by assumption IH-31 as well 
the linearity of the sharp operator, see proposition 12.11 We take the IP-expectation of the 
two right sides. We remark that fa{0, x, uj) must be J^-measurable, so is independent of oj 
and we denote it fa{0, x). In the second equation the integral is split into a Lebesque and 
a stochastic one. Thanks to the assumption IH-4| the stochastic integral is a martingale and 
the expectation is worth zero. Thanks to the Fubini theorem we can exchange the Lebesque 
integral with the expectation. At the end, we find 
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(q,x, 0)fa{0, x) a* 
dA 
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da{fa 



{q, S t , t) S t fa{t, S t , to) 



dt. 



Now, we come back to the equation (|12p . the first term into the expectation is h'{PkL). 
The law of the profit and loss process is worth zero almost surely if we suppose that <; = a. 
Using the assumption IH-51 we can therefore replace h'{PkL) with h'{0). 

We conclude the proof on the relation for the quadratic error r [ E [h {PkL)] ] using the 
independence between the error structure on each parameter and the definition of the 
sharp operator, so we find the relation ([9]). 

The study of the bias is more complicated, we start applying the linearity and the 
closedness of the bias operator: 



(14) A [E [h{PkL)\] = E [A[h{PkL)]\ = E 



h'{PkL) A [P&L] + -h"{PkL) T [PkL] 



where, in the second identity, we have exploited the smoothness of h and the chain rule of 
A. We study the two terms of the previous identity separately. We start studying the term 
depending on the quadratic error. Thanks to relations (I13p and the property of sharp, see 
proposition 12.11 we have 



f dF 

T [PkL] = ^2 0) fa{0, x, u) + 



OA 

da{fa) 



(?, S t , t) fa{t, S t , u) dS t \ T[ai]{ai). 



We remark that the integral into the brackets can be split into a Lebesque and a stochastic 
integral. We can now compute the expectation under P and we find 



E[r [PkL]] = ^2 r [<h] (<h) 



(15) 
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{q, S s , s)fa{s, S s , u)S s ds 



ds 



where we have make use of the properties of the Ito-integrals. We remark in particular that 
E [r [PkL]] ^ T [E [PkL]], since the operator F is bilinear. 
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Now we study the term A[P&zL], we apply the linearity of this operator and we find 



A[PkL] = A [Ffc, x, 0)] + A 



[ Afr S t , t)dS t 

Jo 



A mSr)] 



The last term is worth zero, since the final payoff is completely defined by St and does 
not depend on the volatility estimated by the option seller, see assumptions IH-21 and IH-41 
Thanks to the same assumptions and the closedness of the bias operator A, we have 



A 



f S t , t)dS t ] = [ T A[A(q, S t , t)} dS t . 
Jo J io 



Thanks to the same argument used for the quadratic error, we can take the Gateaux- 
derivatives of F and A with respect to a variation of the volatility along the component fa 
and using the bias chain rule ([2]), we have 



A[F(s, x, 0)] 
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4E 



da(4>i) 
T d 2 A 
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;(q,S t , t) 4> 2 {t, S t , Uj)dS t 



T[ai](ai) 



where we have used that the sum over the index i has a finite number of terms thanks to the 
point 4 of the assumption IH-31 We can take the expectation and, thanks to assumptions 
IH-31 and IH-4| the stochastic integrals are martingales. Therefore, we have 



E [A[PkL\] 



E 



dF 

da{cf>i] 

T 



(16) 
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(q,x, 0) 4>i(0, x) A[oi](ai) 
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ocr{(pi) 

q,x, 0) <p 2 (0, x) r[oi](ai) 



dt A[ai](ai 



o 2 a 



(q,S t , t) (p 2 (t, S t , uj) S t 



dt r[a;](aj). 



We consider therefore the relation (fl"4"|) . Thanks to assumption IH-51 we have h'{P&L) = 
h'(0) and h"(PlkL) = h"(0). We conclude the computation using relations (|15p and ()16p . 
We find the relation (|8|) . The proof ends with the truncated expansion that is a consequence 
of remark |2~T1 and proposition 12, 2| see Bouleau [11] and |12] , 



□ 



3.2 Option Pricing 

In order to interpret this result in finance, we consider that the option seller knows the 
presence of errors in his procedure and wants to neutralize this effect. It is plain that the 
option seller does not control the uncertainties, so the risk related to the space (f2, F , P), 
i.e. to the uncertainty on the parameters, can not be hedged. Furthermore, the relation 
(|10P of the main theorem 13.11 involves that it is not possible to bound this risk, i.e. to 
propose a super- hedging strategy, see Avelaneda and Lyons [26] • Indeed, if we compute 
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the super-hedging price of a contingent claim, we find a very high buy-price and a very 
small sell-price, i.e. the bid-ask spread becomes to large compared with the market spreads. 

This result rises directly from the approach followed. We consider a very large class 
of model and any restrictive hypothesis on the law of the parameters uncertainties. The 
only strong hypothesis is on their magnitude, we have supposed that the uncertainty are 
small compared with the estimated values. We exploit further this property. Thanks to 
that, the probability that the random variable E[/i(P&L)] takes values far from the mean 
is very small and becomes negligible if this distance is big compared with the proper length 
a/ e[h'(0)} 2 ^(a). Therefore, we introduce the following principle to define the price of a 
contingent claim under this uncertainty model. 

Principle P-l (Asset pricing under uncertainty) 

The asset seller fixes a tolerable risk probability a < 0.5 and accepts to sell the option at 
any price F se u, such that 

(17) P {F seU - F + E [P&L] < 0} < a 

where F denotes the cost of the hedging strategy, that is the theoretical price of the option 
without uncertainty, i.e. F = E[&(St)]- 

Before applying this principle to our analysis, we discuss the financial implications of 
this principle. The first remark is that this principle does not give a single price but an 
half- line of possible selling price. Indeed if X is a selling price and Y > X then 

F{Y - F + E[P&zL] <0} <F{X - F + E[PkL] < 0} < a, 

so Y is a selling price too. We call ask price, denoted P a sk 5 the infimum of the set of all 
selling option. We analyze the event {X — F + E [P&L] < 0}. On this event, if the seller 
option has sold the option at the price X, he loses money at maturity, since he has to pay 
F to buy the hedging portfolio and the noise on this hedging strategy conducts him to lose 
money. However, the probability of this event is smaller than a. 

Therefore, the asset pricing principle says that the option seller accepts to sell the option 
at a price X if this price is high enough to guarantee that he loses money with a probability 
smaller than a. It is plain that the parameter a depends on the risk aversion of the option 
seller. However, if he is too risk-adverse, he proposes their contingent claims at a too high 
price, then the buyers can find other traders that offer the same options at a lower price. We 
remark that there is a likeness between this principle and the hypothesis testing in statistics, 
indeed a can be interpreted as the critical probability of the error of the first kind, i.e. "lose 
money in the contract" in our case. The error of the second kind is to propose a too high 
price such that the contract is not signed, i.e. "lose the opportunity to make money". We 
will study the problem of the optimal proposal price in a next paper. We now assume that 
the tolerable risk probability a is fixed and that the option seller chooses the smaller price 
consistent with him tolerable risk probability. 

We also remark that the principle IP-ll defines the purchase price too. As a matter of 
fact, if a trader accept to buy an option, he has to take a negative position on the hedging 
portfolio to cover them, i.e. he has to follow the opposite of the hedging strategy. In accord 
with the principle IP-ll the option buyer accepts to buy the option at any price Pbuy such 
that 

P{-Pbu y + F - E [P&L] < 0} < a, 

that is the relation (|17p where we have changed all signs in order to consider the short 
position, the previous relation can be rewritten as 

P {Pbuy -F + E [P&L] > 0} < a. 

The previous relation and the relation (fT7|) lead us to the following remark: 
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Remark 3.3 (Bid-Ask Spread) 

If all option traders on the market are risk adverse, then they are a difference between 
the best purchase price, denoted Fbid and the best seller price F^k- 

This remark results from the fact that if the trader are risk adverse the tolerable risk 
probability must be smaller that 0.5, then F bid < F + E p [E[P&L]] < F ask . 

The previous principle and the theorem 13.11 have the following immediate consequence: 

Proposition 3.2 (Option prices) 

We assume all hypotheses of theorem \3. 1\ held. We suppose that the option seller follows 
the principle \P-l\ then he accepts to sell the option at any price bigger than 

F ask = F + eA [E[h(PkL)}\ + f [E[h(PkL)]\ M-a 

where M\— a is the (1 — a)-quantile of the reduced normal law, or a function given by Cheby- 
chev's inequality ^ in the conservative case. Likewise, the option buyer accepts to buy the 
option at any price lower than 

F bld = F + eA [K[h(PkL)\] + T \E[h(PkL)]] N a 

Remark 3.4 (Mid price and Bid-Ask spread) 

We remark that the two previous prices are symmetric, since N a -\-M\- a = 0; therefore 
we have 

(18) Fask ± Fbid = F + e A[K[h(PkL)}] 

(19) Fask ~ Fbld = VerpE[/i(P&L)]]A/"i-a 

We emphasize that with our model we can reproduce a bid-ask spread and we can associate 
its width to the trader's risk aversion (the probability a) and the volatility uncertainty (the 
term \/eT \E[h(P&cL)]\). Another interesting point it that the mid-price does not depend 
on the tolerable risk a 

4 Example: log-normal diffusion 

In this section, we give an example of the previous results. In particular, we consider the 
log- normal diffusion, i.e. the Black Scholes model [6]. The underlying follows the SDE 
([5]). In this case, the volatility is a parameter, so we replace them by the same parameter 
multiplied by the identity function. In this model, the assumption IH-31 is simplified, since 
the series expansion of point 4 is replaced by a unique function, i.e. the identity, while 
the points 1 and 3 can be replaced by the more general Hamza condition, see Bouleau and 
Hirsch [TJ. We concentrate our analysis on call option, it is plain that a call option does not 
verify the assumption IH-41 but we have a closed form for all greeks in the case of log-normal 
diffusion and we can check that the Black Scholes pricing formula is C 2 at each time t 
strictly smaller than the maturity T and the vega, i.e. the derivative w.r.t. the volatility, 
vanishes when t goes to T. These properties guarantee that the theorem 13.11 remain true 
even without assumption IH-41 In order to simplify our numerical computation we assume 
that the drift fi = under historical probability P, the computations in the case [i ^ can 
be found on Regis and Scotti |27| . 

In this case, the theorem 13.11 has the following corollary: 

Corollary 4.1 (Bid and Ask prices with log-normal diffusion) 

If the underlying follows the Black Scholes SDE ^ without drift, then the bias and the 
variance of the profit and loss process that hedges a call option verify the following equations: 
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(20) A[C(x, K, T)} 

(21) T[C(x, K, T)] 
where 
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cmd C(x, i^, T) denotes the call of strike K, maturity T and underlying that quotes x. 
Moreover we have the following bid and ask prices: 



(22) C ask (x, K, T) = xNidx) - KN{d 2 ) + V~ex 
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(23) C bid (x,K,T) = xM{d l )-KM{d 2 ) + V~ex 
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Proof: We know that the price of a call on Black Scholes model is given by C(x, K, T) = 
F(o~o, x, 0) = xM(di) — KJ\f(d 2 ), the two previous equation come directly from a computa- 
tion of A\C(x, K, T)\ and T[C(x, K, T)\ using theorem 13,11 for more detail see Scotti |28) . 
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4.1 Analysis of uncertainty impact 

We now analyze the correction on the pricing formula due to the presence of an uncertainty 
on volatility. We have the following corollaries that can be proved with easy computations, 
we analyze only the mid-price, i.e. the average price between bid and ask. We also use the 
notation 



A 



a VT) =2a VT 



?=O"0 



?=O"0 



Corollary 4.2 (Delta and Gamma correction) 

The bias on the call price due to the uncertainty on volatility verifies 



(24) 



dA[C] _ di A[C] 
dK 



e 2 



id 2 

2 a l 



Ka VT 2K UqT V2tt 



(di + d 2 ) r 



Moreover, the bias and its first derivative w.r.t. K are positive (resp. negative) at the money 
if r r S (o'oVTj > ( < )j°~o T. Finally the bias is convex on K at the money if and only if 

?s / - ^+4^+32 
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Corollary 4.3 (Time evolution of the correction at the money) 

The cross derivative of the bias two times with respect to the strike K and one time with 
respect to the maturity, i.e. d ^-^ T is positive if and only if 



r? s (a y/T) > 



1 qgT(qgr-4) 2 + 128 
4 16 + o% T 2 



For the numerical computations see Scotti |28| . The two previous corollaries have an 
interesting consequence. 

Proposition 4.4 (Smile on implied volatility) 

It exists an interval ]a, b[ such that V r^ s (aQ y/T) s]a, b[ the implied volatility computed 
using the price given by corollary \4-l\ is a convex function around the money. Furthermore, 
the second derivative of the implied volatility with respect to strike is a decreasing function 
of the maturity T. That is to say that the implied volatility is smiled and this effect is more 
accentuated for short maturities than long one. 

Proof: We start fixing rf s (o~q\/T) = |<TgT. Thanks to corollary 14. 2| we have that 
A[C] = at the money with its first derivative with respect to the strike. Besides the 
same corollary assures us that the bias A[C] is strictly convex around the money thanks 
the continuity of the derivatives. For the same value of s (ffo VT) we easily check that 
the cross derivative of the bias two times w.r.t. K and one time w.r.t. T is negative at 
the money. We argue that the call price at the money is equal to the Black Scholes price 
whereas around the money the call price in our model is bigger than the Black Scholes 
price. Giving that the Black Scholes formula is used to find the implied, we conclude that 
the implied volatility is strictly convex around the money and that this convexity is more 
marked for short maturities. We now remark that the implied volatility belogns to C°° 
as function of the strike K since it is a composition of infinite differentiable function. We 
conclude that, for each sufficiently small neighborhood of the money, it exists an interval 
]a, b[, including ja 2 T, such that the implied volatility remains strictly convex and with 
convexity decreasing with the maturity. 

□ 
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